Quantum limited velocity readout and quantum feedback cooling of a trapped ion via 

electromagnetically induced transparency 



P. Rabl, V. Steixner and P. Zoller 

Institute for Theoretical Physics, University of Innsbruck, and 
Institute for Quantum Optics and Quantum Information of the Austrian Academy of Science, 

6020 Innsbruck, Austria 



in 
o 
o 

C 

(N 



> 

in 
oo 

vo : 
o . 

in ■ 
o ■ 

^ : 

9 L,: 
^— > ■ 

G ■ 

s : 
cr 



13 



We discuss continuous observation of the momentum of a single atom by employing the high 
velocity sensitivity of the index of refraction in a driven A-system based on electromagnetically 
induced transparency (EIT). In the ideal limit of unit collection efficiency this provides a quantum 
limited measurement with minimal backaction on the atomic motion. A feedback loop, which drives 
the atom with a force proportional to measured signal, provides a cooling mechanism for the atomic 
motion. We derive the master equation which describes the feedback cooling and show that in the 
Lamb-Dicke limit the steady state energies are close to the ground state, limited only by the photon 
collection efficiency. Outside of the Lamb-Dicke regime the predicted temperatures are well below 
the Doppler limit. 



PACS numbers: 42.50.Gy, 3.65.Ta, 42.50.Vk, 42.50.Lc 



I. INTRODUCTION 

Quantum feedback control employs the strategy of act- 
ing on a system based on measurement data obtained 
by continuous observation of the quantum system of in- 
terest, thus achieving control of quantum dynamics and 
preparation of particular quantum states 0, 0, LJ Hi ■ A 
prerequisite of developing quantum feedback control is 
the realization of quantum limited measurements. Quan- 
tum optical systems and, more recently, mesoscopic sys- 
tems have taken a leading role in achieving these require- 
ments towards demonstration of feedback control in the 
laboratory [1 II II II [3 13 13 El 

Motivated by the remarkable experimental progress 
with trapped ions, we will develop in the present pa- 
per a theory of quantum feedback cooling of a single ion, 
based on a continuous readout of the atomic velocity via 
dispersive interactions with laser light. The idea is to 
devise an (in principle) quantum limited measurement of 
the velocity by employing the strong detuning (and thus 
velocity) dependence of the index of refraction of a driven 
A-system near the atomic dark state. These dark states 
are coherent superpositions of two atomic ground states 
which do not couple to the excited atomic state, which 
leads to strong suppression of dissipative light scattering. 
This is the same feature which underlies recent studies 
of electromagnetically induced transparency, slow light 
in atomic gases and quantum memory of light in atomic 
ensembles. The present setup of dispersive readout of the 
atomic velocity complements and is in contrast to ongo- 
ing experiments of quantum feedback cooling of a single 
two- level ion in front of a mirror |l4l ITfij , where the posi- 
tion of the ion is continuously monitored by emission of 
light into the mirror mode, as analyzed theoretically in 
our recent publication |16|. 

Our discussion of quantum feedback cooling of a single 
trapped ion in a strongly driven atomic A-system builds 



on, and connects various well-developed topics in atomic 
physics and quantum optics, as well as continuous mea- 
surement and quantum feedback theory. Thus we find 
it worthwhile to present both a brief review of the back- 
ground material and physical key ideas underlying the 
present work, as well as the main results of the paper in 
Section m In Section ITTT1 we give the technical details of 
our model and derive the equations for the measured sig- 
nal and the conditioned evolution of the atomic motion 
for a weakly excited atom. A Wiseman-Milburn-type pj 
master equation for feedback cooling and the resulting 
temperatures will be discussed in Section llVl Finally, in 
Section[V]we make a connection to EIT-laser cooling and 
describe combination of feedback and laser cooling. The 
details of the adiabatic elimination of the internal atomic 
states are given in the Appendixes and [B] 



II. OVERVIEW AND SUMMARY 



In this section we present an overview of the concepts 
and the main results of this paper. Our emphasis will be 
on explaining the basic physics and strategy behind our 
quantum feedback scheme, and providing references to 
later sections where the mathematically inclined reader 
can find the details of the derivations. In Subsection III Al 
we will briefly review EIT and discuss the dependence 
of the index of refraction on the ion momentum near 
the dark state resonance. Continuous read out of the ion 
momentum using homodyne detection will be formulated 
in Subsection lll Bl The main results of the present paper 
are the equations for feedback cooling in Subsection lll CI 
and the predictions for the final temperatures and cooling 
rates in lllDI in particular also in connection to EIT laser 
cooling of ions [3 13 E3 
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FIG. 1: a) Atomic A-system with two ground states \g) and 
r). The transition \r) — ► |e) is driven by a strong resonant 
light field with Rabi frequency Ol, while the probe field E p 
couples \g) — » |e) with detuning A p . b) Susceptibility x(A P ) 
of a A-system as a function of the detuning of the probe field, 
E p . At the center of the transparency window the real part of 
the susceptibility (solid line) representing the index of refrac- 
tion of the medium exhibits a steep slope, while the imaginary 
part (dotted line) vanishes. 



FIG. 2: Schematic setup for the continuous observation of the 
atomic momentum. The trapped atom is driven by a strong 
laser to create the transparency effect for the probe field, E v . 
Modulations of the probe light are detected in a homodyne 
measurement, i.e., by mixing E p with the strong field of the lo- 
cal oscillator (LO). For feedback cooling, a force proportional 
to the measured signal, I c (t), is applied on the atom. Alter- 
natively to the lens system, the atom can be placed inside a 
running wave cavity, e.g., to enhance the photon collection 
efficiency, e. 



A. Electromagnetically Induced Transparency 

The phenomenon of EIT is related to a quantum in- 
terference effect which in its simplest form can be ob- 
served in a three level atom (for a review see: |2*cl Ell ] 
and references therein). To discuss this effect we con- 
sider an atom with the internal states \g), |e) and |r) 
in a A-configuration as shown in Fig. ^ The transition 
between the states |r) and |e) is driven by a strong, res- 
onant laser, while a second light field, the probe field, 
couples the ground state to the excited state. We denote 



susceptibility has the characteristic behavior [2l| 



by A p 



the detuning of the probe field from 



the atomic resonance. The atomic Hamiltonian is then 
given by 



H A = HA p \g)(g\ + -(n L \e)(r\+g\e)(g\ + h.c.) , (1) 

where l and g are the Rabi frequencies of the laser and 
the probe field respectively. At the two photon resonance, 
Ap = 0, the Hamiltonian, H A has an adiabatic eigenstate 



with zero energy, a so-called "dark state" , 



\D) 



tt L \g)-g\r) 



(2) 



For an atom in the state \D), the excitation from the 
states \g) and |r) destructively interfere and the atom 
decouples from the light. We note that an atom in a 
dark state involves no excited state population, and is 
thus immune to decay from the excited state. 

The existence of such a dark state leads to remarkable 
properties of the index of refraction. For weak probe 
fields, the propagation can be discussed in terms of the 
linear susceptibility, x( Ll 'p)- The real part of xi^p) IS 
related to the refractive index by n — 1 + Re(x(w P ))/2, 
while the imaginary part, Im(x(wp)), is proportional to 
the absorption coefficient. For an ensemble of three level 
atoms where both, \g) and |r) are long-lived states, the 



iA T 



\{l L \ 2 /4 + iA p (T + iA t 



(3) 



where T is the decay rate of the excited state. Fig. ^ 
shows the dependence of x 011 the detuning, A p . Around 
the dark resonance A p = we have a steep slope of the 
refractive index (solid line) which leads to a slow group 
velocity of the probe field ("slow light") while absorp- 
tion is strongly suppressed (dotted line), giving rise to 
"electromagnetically induced transparency" . The width 
of the transparency window as well as the variation of 
the refractive index depend on and can be controlled 
by the laser field. We note the suppression of absorption 
at the dark state resonance (dotted line in Fig.^Jj). 

Consider now a single ion in a A-configuration (Fig.[3i) 
moving in a trapping potential. We will only consider 
the ID motion along the propagation direction of both 
the probe and dressing laser beams. If we adopt for the 
moment a classical description of the ion motion with 
z(t) the ion trajectory, then the internal dynamics of the 
ion can again be described by the Hamiltonian JQ) with 
a Doppler-shifted probe detuning 



Ap(i) 



«(*): 



where k p and kh the wave vectors of the running probe 
field and the dressing laser field, respectively, and v(t) = 
z(t) is the ion velocity. Thus for a resonant probe field, 
Ap = 0, the change in the index of refraction is a linear 
function of the atomic velocity or momentum with the 
steep slope given in Eq. 

This suggests the strategy to measure the momentum 
of the ion continuously by monitoring the phase shift 
due to the varying index of refraction. We note that - 
while the index of refraction of a single particle is small 
- EIT will strongly amplify the sensitivity to the veloc- 
ity. At the same time dissipation due to light scattering 
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is strongly suppressed within the transparency window. 
These arguments can be easily adapted to a situation 
where the atomic motion is quantized. 

B. Continuous observation of the momentum of 
the trapped ion 

We turn now to a formulation of the continuous read 
out of the ion momentum as outlined in Fig.|3 The idea 
is to measure the momentum of the ion via the phase 
changes of the probe beam as described in the previous 
subsection. The phase of the probe beam can be deter- 
mined by homodyning, i.e. by mixing the probe beam 
with a local oscillator and measuring the homodyne cur- 
rent, I c (t). 

The state of the observed system (the moving ion) is 
described by a conditional density operator /i c (i), which 
represents the observer's knowledge of the current state 
of the system for a given record of the measured signal, 
I c (t). We will show in Sec. HHI that after adiabatic elimi- 
nation of the excited state of the (weakly driven) ion the 
measured homodyne current has the form, 

J c (t) = 2eT {p) c + v^£(t) • (4) 

The current is the sum of two terms. The first contribu- 
tion shows a linear dependence on the conditional expec- 
tation value of the momentum operator (p) c = Tr{p, c p}. 
Thus by measuring I c (t) we learn the momentum of the 
moving ion. The second contribution describes a shot 
noise term with £(i) a white noise Gaussian process. The 
signal strength is determined by the rate Tq. It is related 
to the slope of the refractive index x & t A p = 0: an ex- 
plicit expression is given in Sec. IIIII Ea. i|39[l below. The 
parameter < e < 1 takes into account the collection ef- 
ficiency of the scattered photons (where in the ideal case 
e = 1). Eq. 10} is derived under the assumption that 
is large compared to typical Doppler detuning, Ad, and 
is valid on time scale which is slow compared to fi^ 1 . The 
signal is maximized for the local oscillator phase, = 0. 

According to continuous measurement theory applied 
to homodyne detection, the conditional density operator 
/x c (t) is updated upon observation of the current I c (t) 
following the Ito equation, 

dficif) = — iv[a*a, fj, c (t)]dt (5) 
+C M ^c(t) dt + T/eToH[p\n c {t) dW(t) . 

This equation will be derived in Sec. IIHI Eq. Ij38(l . 

The first term in this equation describes the free evo- 
lution in the ID harmonic trap where v denotes the trap 
frequency, and a {a)) the destruction (creation) opera- 
tors, respectively. The effects of the continuous observa- 
tion appear in the second and third term of Eq. © . 

The superoperator Cm determines the back action of 
the measurement setup on the atomic motion. In the 
Lamb-Dicke limit r\ = 27rao/A p <C 1, where the extension 
of the atomic wavepacket (size of the harmonic oscillator 



ground state do) is much smaller than the wavelength of 
the light, X p , it has the form, 

p 

The action of Cm tends to diagonalize the density op- 
erator in the eigenbasis of the (measured) operator p. 
By comparing the decoherence rate, To, with the signal 
strength, eFo, we see that for e < 1 the measurement 
is not quantum limited, i.e., more noise is added than 
required by quantum mechanics . Although the mea- 
surement does not reach the quantum limit, the back 
action is still minimal for a given collection efficiency. 
In the third term of Eq. JSJ we introduced the notation, 

Ti,[c]fx = c/i + fxc — (c+ c)/i . 

This term describes the observer's knowledge of the cur- 
rent state of the system and therefore depends on the 
measured signal. The Wiener increment dW(t) is for- 
mally related to the signal noise by dW(t) = £(t)dt. 

In summary, Eq. (0J for the homodyne current I c {t) 
and the evolution equation (J5J for the conditional den- 
sity matrix constitute the basic equations of continuous 
observation of the momentum of the ion via homodyne 
detection. 



C. Quantum Feedback Cooling 

The information on the atomic momentum contained 
in the signal I c (t) I0J can be used to act back on the sys- 
tem. Here we are interested in cooling the atomic motion 
by u sing the feedback strategy known as "cold damp- 
ing" [23( . The idea is to apply a force on the atom which 
is proportional but opposite to its momentum. This force 
creates an effective friction for the atomic motion and, 
therefore, leads to a dissipation of kinetic energy. 

In our setup the measured signal is already propor- 
tional to the average momentum, (p) c , and can be am- 
plified and fed back directly. Thus we consider a feedback 
Hamiltonian of the form, 

H fb (t) = ^I c (t-r)z, (7) 

where G is the dimensionless gain factor, and r is the 
finite delay in the feedback loop. Note that r > 0, so 
that Hfi, acts after the measurement. Eq. (jSJ) with the 
feedback Hamiltonian added provides us with a feedback 
equation describing the time evolution of the system. The 
goal is now to average this equation over the Gaussian 
white noise 

A general theory for direct quantum feedback has been 
first discussed in a seminal paper by Wiseman and Mil- 
burn Q]. In particular, they have shown how to average 
the quantum feedback equation in the limit r — > + . In 
our case this assumption implies that the time delay of 
the feedback is small on the scale of the (adiabatically 



4 



^ 6 
4 

2 



0.2 0.4 



0.6 

G 



0.8 



FIG. 3: The figure shows the dependence of the steady state 
energy as a function of the feedback gain, G. The results are 
calculated in the Lamb-Dicke limit and for F = 0.01 v. The 
three curves are plotted for the parameters, e = 0.1 (solid 
line), e = 0.05 (dashed line) and e = 0.01 (dotted line). 



eliminated) system evolution, a condition which is real- 
istic in the present context. Adopting this formalism we 
will derive in Sec. HVI a master equation for the uncondi- 
tioned density operator, fj,(t) = E\p c (t)]. We obtain (for 
A p = A L = 0) 

fi = — iv[a'a, fi\ + £m(i 

t G r s, .u -i r G Y r 11 (8) 

Z 06 

The first line of this equation describes the free evolu- 
tion and the measurement back action, Cm Eq. ®- The 
terms in the second line of Eq. © include the effects 
of the feedback loop. While the term proportional to G 
causes the expected damping of the motion, the second 
term leads to a diffusion of the atomic momentum. This 
diffusion originates from the noise in the measured cur- 
rent which is also amplified and fed back to the system. 

The generic dependence of the steady state energy on 
the feed back gain calculated from Eq. (JSJ is plotted in 
Fig- El The curves show the expected minimum as a func- 
tion of G, at the point where the noise of the feedback 
loop starts to dominate over the damping force. More 
detailed results will be presented in the following subsec- 
tion. 



D. Results: Quantum Feedback vs. EIT Laser 
Cooling 

In general, the interaction of atoms with light always 
leads to some form of laser cooling or heating. In the 
resonant case A p = At = (discussed above), which is 
required to measure the atomic momentum, heating and 
cooling rates are equal and cause the diffusion described 
by Cm- 

By detuning the lasers away from the resonance, At ^ 
0, the atomic susceptibility, x(A p ), and therefore the ab- 
sorption properties become quite asymmetric (see Fig.0J. 
This asymmetry is exploited in EIT laser cooling (ELC), 




FIG. 4: Linear susceptibility, \(Ap) 
A L /T = n L /T = l. 



A P /r 



in arbitrary units for 



where the absorption on the red sideband (a phonon is 
removed from the motion) is much more likely than on 
the blue sideband (a phonon is added to the motion). A 
detailed discussion of ELC can be found in Ref. 0] . 

In Sec.[V]we show that in the Lamb-Dicke limit we can 
derive the extension of Eq. (JHJ) which includes quantum 
feedback cooling as well as ELC. In rotating frame with 
respect to the trap frequency, v, it can be written in the 
form, 



(A_ + A^)T>[a]n + (A+ + A^)V[a f 



(9) 



with 



a^ia' — —a'a^J,— —[ia' a . 



The total cooling and the total heating rate are divided 
into a contribution from the atom-laser interaction, A±, 
and a contribution from the feedback force, A f ±. They 
are given by, 



A+ = 



-Re[/(±i/)] , 



4 6 = r (G^im[r(±^] + ^ 



The relation between the four rates is determined by the 
function I(y) which is defined as 



I{ V ) 



9. 4 



2Yu 2 {ti 2 - Av(v - A L )) + i2Tv 



(10) 



Note that the rates A± are proportional to the imaginary 
part of the atomic susceptibility at the sideband frequen- 
cies, x( w p T v )i ( see Fig- EJ- For a general detuning At, 
the measured signal is no longer proportional to (p), and 
therefore the phase of the local oscillator, <j), appears in 
the feedback rates. 

In the basis of harmonic oscillator states the master 
equation © has the form of a standard rate equation for 
the trap occupations (p n = (n\p\n)) familiar from laser 
cooling, 



Pn =(A_+ A f _ b )[(n + l) Pn+1 



np r , 



+ (A + + A f + b ) [(n - l)p B _ x - (n + l)p n ] 



(11) 



5 



b£: 

o 




FIG. 5: Steady state occupation number, n — {a/ a), as a 
function of the detuning, Al- The solid line shows the result 
for optimized (j> and optimized gain, G, while the dashed line 
shows the result without the feedback loop. The curves are 
plotted for the parameters: e = 0.05, Ql/T = 0.8 and v/T = 
0.1. 



which predicts in steady state a Bose-Einstein distribu- 
tion with a mean occupation number 



-4, 



(A_ - A+) + (A f _ b - A{ b ) 



(12) 



We now turn to the discussion of results for the case of 
pure feedback cooling, and combined feedback and laser 
cooling: 

Pure feedback cooling. We first reproduce the results 
of Eq. © by setting A p = A L = 0. Then A± = T /2 
and the cooling of the atom is attributed to the feedback 
mechanism. The minimal energy is reached for G = \fAe 
and 4> — 0. It is given by 



hp 
T 



(13) 



This expression shows that the final temperature is only 
limited by the collection efficiency e. In the theoretical 
limit, e — ► 1, it approaches the ground state energy, Kv/2. 

Feedback cooling and ELC. When we tune away from 
the resonance, A p = ^ 0, we obtain a difference 
in the laser cooling rates, A^ ^ A + . In Fig. [S] we com- 
pare the final temperatures for the optimal feedback gain 
with the case of pure ELC. For blue detuning, A^ > 
the mechanism of ELC sets in and cools the atom close 
to the ground state. Although the addition of the feed- 
back loop always leads to even lower temperatures, its 
effect can be neglected because for the present parame- 
ters ELC already provides efficient ground state cooling. 
For red detuning A^ < the absorption spectrum (see 
Fig. 2J) is reversed and the atom is actively heated by 
the light absorption. In this case a steady state is only 
reached when the feedback cooling dominates over the 
laser induced heating. 



III. THE MODEL 

In this section we present a detailed description of 
our model for the measurement setup which is shown 
in Fig. |21 The system of interest is the three level atom 
which is confined by an external trapping potential. This 
atom is illuminated by a strong laser field to create the 
transparency effect for the probe field. The outgoing 
probe light is mixed with a strong local oscillator to per- 
form a homodyne measurement to detect linear shifts of 
the field. 

To describe the dynamics of the atom as well as the 
detection of scattered field, we start with the total Hamil- 
tonian, 



H = H A + H A - 



EM 



H 



EM 



(14) 



It is the sum of the Hamiltonian for the external and 
internal states of the atom, Ha, the free Hamiltonian for 
the electromagnetic environment, Hem, and the coupling 
between the atom and the electromagnetic field, Ha-e- 
For the internal level structure we consider a A- 
configuration as shown in Fig. A classical laser field 
with frequency lul drives the transition between the ex- 
cited state, |e) and the second ground or metastable 
state, |r). The Rabi frequency for this coupling is de- 
noted by For the external dynamics of the atom 
we restrict ourselves to a one dimensional model, i.e., 
we assume that the atom is strongly confined in the x 
and y directions. Along the z-axis, which coincides with 
the propagation direction of the probe beam, the atom 
is trapped by the external potential, V(z). Although it 
is not essential for the following discussion, we further 
assume that V(z) is harmonic, with a trap frequency, v. 
This assumption allows us to introduce the dimcnsionlcss 
position and momentum operators, z := (a + a* )/v / 2 and 
p := z(<jt — a)/V2, where a and denote the usual an- 
nihilation and creation operators. With the definition of 
the external Hamiltonian, He = hva)a, and the notation 
o~ij = \i)(j\ the atom Hamiltonian is then given by 



H A = H E 



hio e \e)(e 



■ hw r \r) (r\ 



h.c) 



(15) 



Here aj e ,r denote the eigenfrequencies of the correspond- 
ing states |e) and \r). The Lamb-Dicke parameter is de- 
fined as, r/ r = k^y/h/mv, where m is the mass of the 
atom and k^ is the wave vector (projected on the z-axis) 
of the laser field. 

The electromagnetic environment consist of a three di- 
mensional set of plane wave modes, labelled by their wave 
vector, k and their polarization, A. In terms of the corre- 
sponding annihilation and creation operators, b\{k) and 
b\(k), the free evolution is determined by the Hamilto- 
nian, 



Hem = 



A=l,2 



d 3 khLUj:bl(k)bx{k). 



(16) 
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The electric field of the environment interacts with the 
internal states of the atom via a dipole coupling. Under 
the rotating wave approximation the interaction Hamil- 
tonian is 

i?A-B = -fiegE + (z)a eg - p er E + (z)a er + h.c. , (17) 
with the standard expression for the electric field, 



E+(x) = i / d3k 



A=l,2 



hwi, 



2e (27r) 



■e x (k)b x (k)e lkS . (18) 



In an experiment the lens system defines a certain spatial 
mode function for the probe beam. To describe the ho- 
modyne detection of the probe field, we divide the total 
electric field into the field of this particular mode, E p , 
and a remaining set of modes, orthogonal to E p , 

E+(x) = £+{$) (19) 

In our model we approximate E p by the one dimensional 
field, 



E+(z)= l £ p 



^(kpZ—ujpt) 



oo 

ie p I dka(k) e ikz b p (k) . 
o 



The first part in this expression describes the coherent 
field of the incoming probe beam. Note that with this 
definition of the operator, E pi the initial state of the 
electromagnetic environment is the vacuum state. The 
function, a(k), determines the coupling of the atom to 
the dense set of modes, b p (k). It must be adjusted to 
reproduce the correct results of the real mode function. 

The interaction of the atom with the coherent part 
of the probe field leads to transitions between \g) and 
|e), characterized by the Rabi frequency, g — 2\(2 eg £ p \/h. 
In the following include this term in the atomic evo- 
lution and define a new system Hamiltonian, H$ = 
Ha + Ha-em\co\\- In a frame rotating with the laser 
frequencies, oj p and ll>l, this Hamiltonian is then given 
by 

H S = H E - HA p \e){e\ - h(A p - A L )\r)(r\ 

(21) 



if 



, h 9 , 

+ Y ,f 



-'"7s 



'(Tge) ■ 



Here we introduced the detunings A p — lu p — to e and 
Al = wi — (uj e — w r ), and the Lamb-Dicke parameter of 
the probe field, r\ g — k p \Jhjmv. 

The coherent evolution of the driven atom is now de- 
termined by the system Hamiltonian Hg. The coupling 
between the atom and the dense set of modes of E± and 
the non-classical part of E p has two effects. First, it leads 
to an incoherent dynamic of the atomic state. This in- 
cludes the decay of the excited state population as well 
as a diffusion of the atomic momentum due to the ran- 
dom recoil kicks of the emitted photons. Second, the 
coupling changes the state of the electromagnetic envi- 
ronment which is (partially) detectable by the observer. 



A. Master equation 

We first look at the incoherent dynamics of the atom 
and ignore the evolution of the electromagnetic bath. By 
applying the standard Born-Markov approximation and 
tracing over the bath degrees of freedom we obtain a 
master equation for the system density matrix, p. It can 
be written in the standard form, p — Cp, with a Liouville 
operator EHHm, 

Cp = --[Heffp - pHl s ] + Jg{(JgeP<Jeg) + Jr^reP^er) ■ 

(22) 

The effective Hamiltonian, H c s = H$ — ihT/2\e) (e|, in- 
cludes the unitary evolution of the atom as well as the 
decay of the excited state population with a total rate, 
r. The two "recycling" terms are defined by 

Jj=r,g(p) = T J J duN(u) e -^ iu pe l ^ u . (23) 

The rates T g and T r denote the decay rates into the cor- 
responding states, \g) and |r). The dipole distribution, 
N(u) — |(l + u 2 ), with u = cos(ip), determines the prob- 
ability for emitting a photon under a certain angle, if, 
with respect to the z-axis. 

Master equation l|22|) describes the full dynamics of a 
driven atom in a A-configuration. The external and in- 
ternal degrees of freedom are coupled via the position 
dependent interaction with the electromagnetic field. In 
general, the recoil kicks of the emitted photons described 
by Jg, r lead to momentum diffusion and a heating of the 
atomic motion. For an appropriate choice of laser detun- 
ings the heating can be compensated by photon absorp- 
tions (laser cooling). In this paper, we follow a different 
approach where cooling is provided by an external feed- 
back force. 



B. Continuous homodyne detection 

In a next step we describe the homodyne detection of 
the probe field, E p . Here we follow the standard theory of 
homodyne detection (see, e.g. Ref. [2(| |2?], 13 ) to derive 
the relevant equations for the model specified above. 

After the interaction with the atom, the outgoing 
probe field is desribed by the Heisenberg operator 



(t). Using the input-output formalism |26| it is re- 



lated to the incoming field, b Pt i n (t), by 



^p.out 



(t) 



(24) 



where c D 



-iri p z 



a ge denotes the atomic "jump opera- 



tor" which couples to the probe field. To relate our model 
to the real experimental setup, we set 7 = eT g , where the 
collection efficiency e determines the fraction of photons 
which are scattered into the mode of the probe field. 

In homodyne detection the outgoing field (|24|l is mixed 
with a strong coherent field of the local oscillator. When 
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the transmittance of the beam splitter is close to one the 
field operator at the position of the detector is 

bd(t) = y/iPeVe-*"'* + b p . in (t) + ^c p (t) . (25) 

Here, j3 and <f> denote the real amplitude and the phase 
of the reflected part of the local oscillator. Note that 
the total probe field defined in Eq. (|2U[) is the sum of 
a classical and a quantized contribution. The classical 
part of bp^ n {t) can simply be absorbed in a redefinition 
of (3e 1 ^ . The operator for the homodyne current is, 

h(t)= lim (b\(t)b d {t)- 1 p 2 ) /p. (26) 

The measured signal I c (t) is then defined as the outcome 
of the continuous measurement of the current operator, 
Ih(t). Using the results form the theory of homodyne 
detection |26j . we obtain 

I c {t) = eTgidpe-^ + cle^) c (t) + • ( 27 ) 

The unconditioned evolution given by master equa- 
tion ilL'21) and the measurement record, Ic{t), determine 
the evolution of the conditioned density operator, p c (t). 
Following Ref. [H we obtain, 

d Pc (t) = Lp c (t)dt + ^T g H[c p e- l % c {t) dW(t) . (28) 

Eqs. 1)27(1 and (|28[) represent a full description of the 
conditioned dynamics of the atom under continuous ob- 
servation and serve as the starting point for the following 
discussion. 

C. Adiabatic Elimination 

The current I c (t) as given in Eq. (|27|l is still a function 
of the coupled external and internal states of the atom. In 
the following we show that for a weakly excited atom we 
can eliminate the dynamics of the internal states, and the 
measured signal becomes a linear function of the atomic 
momentum as given in Eq. |0J. In addition we derive the 
resulting back action of the measurement on the motional 
state of the atom. 

As already noted in Section [nj the phase shift of the 
probe light is a linear function of the atomic momentum 
as long as the typical Doppler detuning, Ad = iyrj p (p), 
is small compared to the width of the transparency win- 
dow, Ql. Therefore, we can apply perturbation theory 
in the parameter Ad/^l to derive an effective equation 
for the external state. As a first step in our calculation 
we perform a unitary transformation, 

jj = e «?p2| e )<e| e -«7S|r)(r| ; ( 2 9) 

where we set fj := r/ p — t)l. In the new basis the system 
Hamiltonian H$ = U'HsU is given by 

H s = H E -h (A p - vi lp p - vr$/2) |e)(e| 

-h(A p -A L + vfjp-vf/2)\r){r\ (30) 
MIl , , fig i x 

H 7T~ \ a er + O re ) + — {(J eg + (J ge ) . 



As in the classical case (see Section Hi A|) the position de- 
pendence of the atom-laser coupling is transformed into 
a frequency shift for the states |e) and |r). In addition 
to the Doppler detunings, vrj p p and vfjp, the internal 
states are also shifted by the appropriate recoil frequen- 
cies. They account for the fact, that each absorption or 
emission of a photon also transfers kinetic energy to the 
atom. 

For the particular choice of laser detunings, A p = 
+ vfj 2 /2, and in the absence of a trapping potential 
the Hamiltonian Hg has a dark eigenstate, — \p — 
Q,D). Here \p = 0) is the zero momentum eigenstate and 
\D) denotes the internal dark state, 

\D) = (g\r)-n L \g))/n, (31) 

with il = \/Q 2 L + g 2 - In the following we assume that 
this relation between the detunings is fulfilled. The sys- 
tem Hamiltonian can then be written as 

Hs^He + H!+ H int , (32) 

such that He and Hi act on the external or the internal 
states only, while Hi nt describes the coupling between 
them, 

H iat = hvr) p p\e){e\ - hvfjp\r)(r\ . (33) 

With the definition A = A p — vr] 2 /2 the Hamiltonian of 
the internal states, -H/, reduces to the one of a driven 
A-system at the two photon resonance, 

TiVIl hg 
Hi = -nA|e)(e| H — (a er + a re ) + — (a eg + a ge ) . 

(34) 

In the limit where the external and internal degrees of 
freedom decouple, H int — > 0, the conditioned dynamics 
determined by Eq. I|28|l leads to a relaxation of the atom 
into the state 

p c (t^oc)=fi ®\D)(D\, (35) 

where po is an undetermined state of the external de- 
grees of freedom. The interaction with the atomic mo- 
tion, H m t, or to be more precise the term hufjp\r)(r\ cou- 
ples the state |D) to the bright (internal) eigenstates of 
Hj, |+) and |— ). They are given by 

|+) = cos(9)\e)+sm(9)(g\g)+n L \r))/n, 
|-) = sm(d)\e)-cos(9)(g\g)+n L \r))/n, 

where the mixing angle 8 is defined by the relation 
tan(6>) = rt/(\/fl 2 + A 2 — A). Theses two states are sep- 
arated from the dark state by the energies 

hn ± = -h(A T Vtt 2 + A 2 ) /2 . (36) 

Therefore, for small Doppler shifts the population in the 
bright states is of the order of A 2 D / ti± . 

In the following we consider the limit where the eigen 
frequencies of the bright states, f2±, are much larger than 
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the typical Doppler detuning, Ad, as well as the fre- 
quency of the trap, v. The first assumption says that 
the atom is only weakly excited and the total density 
operator is well approximated by, 



p c (t)~Po(t)(8|I>)(D| 



(37) 



The second condition, Ql 3> v, ensures that the internal 
state of the atom adiabatically follows the evolution of 
the atomic momentum. If both conditions are satisfied 
we can adiabatically eliminate the population in the the 
bright states and derive an effective equation for the evo- 
lution of the motional state, p,o(t). The details of this 
calculation are summarized in Appendix ^ Finally, we 
revert the unitary transformation, U l|29|l , and trace over 
the internal states. For the resulting conditioned density 
operator of the motional state, fi c := Tii{Up c U^}, we 
obtain the stochastic master equation, 



dfj, c (t) = — i[va)a — A 2 Ap 2 , fi c (t)]dt + £M/i c (t) dt 
+ ^eVT g H[pe~ l % c {t)dW(t) , 

and the homodyne current 



I c (t) = 2e\ 2 T g cos(0)(p) c + yje\ 2 r g £(t) . 



(38) 



(39) 



In these two equations we defined the parameter, A = 
2gfji//Q 2 , and the average is take with respect to the 
conditioned motional state, (•)<. = Tr{/i c -}. The mea- 
surement back action has the form 



To 



-1 Ur 



(40) 



of the "cold damping" feedback strategy. By applying the 
theory of direct quantum feedback [IJ we derive a master 
equation for the unconditioned state, fi(t) = E[/j, c (t)]. 

For the feedback cooling we consider the measurement 
setup as described in the previous section. In addition, 
we apply a force on the atom which is proportional but 
opposite to the measured current. For a single trapped 
ion, such a feedback loop can be realized by converting 
the homodyne current into a voltage difference between 
two trap electrodes. The effect of the feedback loop on 
the system evolution can be written in the general form, 



fic(t)\fb = Ic(t - T)fCfl c (t) . 



(42) 



The time delay of the feedback loop, r, can usually be ne- 
glected compared to the timescale of the atomic motion, 
v . Nevertheless, for the derivation of the final master 
equation a finite value of r is important to obtain the 
correct operator ordering Q. To implement the idea of 
"cold damping" , we consider the feedback superoperator, 



Kfx = -i— [z,fi\ 



(43) 



where G denotes the dimensionless gain factor. Note that 
with this definition of IC, the frequency scale of the feed- 
back contribution is again of order To. 

It has been shown in Rcf. that Eq. O must be 
interpreted in the Stratonovich sense. To be compatible 
with Eq. (J5J we convert it into the Ito-type equation, 



with r n = A 2 r and 



dfx c (t)\ fb =r (2e(p) c (t - t)JC + |/C 2 ) fi c (t)dt 
+ ^/7r~ Kii c (t)dW{t - t) . 



(44) 



r 



dttJV , (u)e -<1 w 1 < 1i - 1 W w * ( "~ 1) 



(41) 

Discussion. The results given in the Eqs. I|38|l and 139|) 
are valid in the limit of a weak probe field, g <C fiz,- In 
that case the dark state, \D), almost coincides with the 
ground state, \g), and the signal strength is maximized 
for a given strength of the measurement back action, To, 
(see Appendix 0. The signal, J c (i), can further be opti- 
mized setting cf> = and by choosing atomic states with 
a small branching ratio r r /r s . The remaining difference 
between A 2 F g and Tq can be absorbed into the definition 
of e. For A = we then end up with Eq. Q and Eq. © 
as given in Section ITTbI 



In this form we already see that the noise in the current 
I c (t) leads to the diffusion term K? . We can now add 
Eq. (|44|l to the conditioned evolution given in Eq. J5J) and 
obtain the full stochastic dynamics of the atom under the 
action of the feedback loop. 

To derive a master equation which is independent of 
the measurement outcome, we perform an ensemble aver- 
age over the stochastic process, £(t) and obtain the evo- 
lution of the unconditioned density operator, /i = E\p c ]. 
When taking the average we must keep in mind that al- 
though E[dW{t)] = 0, the Ito increment, dW(t — t), is 
not independent of fJ, c (t), e.g., E[dW(t — r)/x c (f)] ^ 0. 
The way to perform the average in the limit, r —* + , 
can be found in Ref. [2^| . By following this procedure we 
end up with the master equation JHJ given in Section ITT1 



IV. FEEDBACK COOLING 

As already mentioned in Section[n]the goal of the con- 
tinuous momentum observation is to use the information 
in the signal to manipulate the motion of the atom, e.g., 
to cool it. In this section we discuss the implementation 



A. Feedback cooling in the Lamb-Dicke limit 

We first look at the solution of Eq. © in the Lamb- 
Dicke limit, T) p ,r)L <C 1. In this limit the recoil kicks of 
the emitted photons can be neglected and the backaction 
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of the measurement <|40ll simplifies to 



(45) 



For a harmonic trapping potential, He = Hisa^a, the 
feedback master equation JSJl is then quadratic in the 
position and the momentum operators. Therefore, the 
final state is Gaussian and we obtain analytic expressions 
for the variances of z and p. The resulting steady state 
energy is given by 



E = — 



hv_ ( GTq 
2 i 2^" 



1 

G 



G 



(46) 



The first contribution in the brackets originates from the 
enhanced uncertainty of the position coordinate as a re- 
sult of the measurement of the momentum operator. If 
the measurement strength, To, is much smaller than the 
trap frequency, v, position and momentum coordinates 
are mixed sufficiently fast and this contribution disap- 
pears. In this limit the optimal feedback gain is given 
by G = \fAe and we obtain the minimal energy given in 
Eq. O. 

In Section El we extend the discussion of the feedback 
cooling in the Lamb-Dicke limit to arbitrary detunings, 
A. Then laser cooling effects play an important role for 
the final temperatures. 



and we obtain an equation for the mean occupation num- 
ber 



(h) = -T \G - D](n) + ^ 



r 2 

— -G + 1 + D 
4e 



(48) 



In this expression, the parameter, D, describes the heat- 
ing induced by the recoil kicks from the emitted photons. 
It is given by 



D = if Q a 



V — V 



■{rfca + r] r rjg) 



r,. 



(49) 

where a = h J N(u)(u — l) 2 du — 7/10. By choosing an 
appropriate atomic level configuration, T g 3> T r and/or 
jjj > ^, its value is only limited by D w fjgOL- In this 
case, and for optimized gain the minimal steady state 
energy is 



En 



j] 2 g a + y^e + j]\h 2 ' 
4e 



(50) 



This expression shows that the minimal energy changes 
from the Lamb-Dicke to the non-Lamb-Dicke regime at 
the parameter values, 4e f=a In the non-Lamb-Dicke 
regime, i.e., for weak trapping potential the minimal en- 
ergy approaches the value, 



B. Feedback cooling beyond the Lamb-Dicke limit 

When the trapping potential is weak, the extension of 
the atomic wavepacket can be of the order of the wave- 
length of the emitted photons. In that case, the energy 
spacing in the trap, Hp, is comparable to the recoil energy, 
En, and recoil kicks from the emitted photons lead to an 
additional diffusion of the atomic momentum. Therefore, 
we must take into account the full expression for the back 
action term Cm- By expanding Eq. (|40|l in the parameter 
r r /r we can write it as 



E-rr\ \ v 



E T 



(51) 



r r ° 



2J„ 



yrc=0 



J r n (P M P)-P 2 M-MP 2 • (47) 



The zeroth order term in the sum corresponds to the 
physical picture where the atom is excited and simply 
decays back to the ground state, \g) w \D). Processes 
where the atom first decays into the state |r), is then 
reexcited again are taken into account by including higher 
order terms in this sum. 

In general, the full expression of Cm leads to a hierar- 
chy of coupled equations for the moments of p which docs 
not break off as in the Lamb-Dicke limit. In the following 
we restrict our discussion to a finite trapping potential 
and consider the limit, Tq <^ v. As mentioned above, this 
ensures a mixing of position and momentum coordinates 
and therefore, an equal reduction of both variances. In 
this regime non-energy conserving terms can be neglected 



Therefore, the limit for feedback cooling is set by the 
recoil energy, Er, divided by e, and temperatures well 
below the Doppler limit, fesTo = hT/2, can be reached. 



V. FEEDBACK VS. EIT LASER COOLING 

In the previous section we focused on the laser detun- 
ings A p ~ A j, ~ with the goal to measure the atomic 
momentum to achieve quantum feedback cooling. As al- 
ready discussed in Section lll Dl in a A-systems EIT laser 
cooling provides an effective mechanism to cool atoms 
essentially to the ground state without any further ex- 
ternal manipulation. In this section we derive a master 
equation which describes both effects, feedback cooling 
and ELC, and discuss the cross-over from pure feedback 
cooling to ELC. 

For the adiabatic elimination of the excited states in 
Section llll Cl and therefore for the validity of the feedback 
master equation © we required, v <C fi±- This assump- 
tion excludes the parameter regime, where ELC achieves 
the lowest temperatures, v ~ Q+ 18]. In this Sec- 
tion we restrict the discussion to the Lamb-Dicke limit, 
VgiVr ^ 1- This allows us to derive a master equation for 
the motional state for arbitrary choice of the parameters 
A, VLl and v. 

We start with the full model for the three level atom 
coupled to the radiation field as introduced in Section lTTTl 
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To optimize the feedback cooling effect and to simplify 
the following discussion we make the assumptions, T g — 
r, and as in the previous sections, g <C f^. Under the 
two photon resonance condition, A p = Aj, = A, the 
system Hamiltonian, Hs, given in Eq. I|21|) can be written 
as 



Hs — He + Hi + H„ 



(52) 



The Hamiltonian H v describes the coupling between ex- 
ternal and internal degrees of freedom. Up to first order 
in the Lamb-Dicke parameters it is given by 



H n 



ihz 



Ml, 



<Tre) + Vg^eg 



<Jge) 



(53) 



The conditioned dynamics of the full atomic density 
operator, p c {t), is determined by the stochastic master 
equation (|28p. As in Section llll Cl the goal is to eliminate 
the internal states and to derive an effective equation for 
conditioned motional density operator, fx c {t). The prin- 
cipal strategy is the same: For vanishing Lamb-Dicke pa- 
rameters the decay of the bright states relaxes the atom 
into the state, p c {t) = p, c (t) <8 \D)(D\. The dynamics of 
Pcif) can be derived by including the coupling Hamilto- 
nian H n in second order perturbation theory. In contrast 
to Section llll Cl wc impose no restrictions on the energies 
of the bright states, f2±, which leads to resonant transi- 
tions for |f2±| w v. Therefore, to guarantee the validity 
of the perturbation theory we require that {H v ) « fjg is 
much smaller than the decay rates of the bright states, 
r+ - T cos 2 (<9) and L_ - Lsin 2 (6>). 

In Appendix[B]we use the stochastic Schrodinger equa- 
tion formalism for the adiabatic elimination of the inter- 
nal states. As a result we obtain the conditioned master 
equation, 

dfic = — i{y + 8) [d^d, p. c ]dt 

+ A-V[a]fj. c dt + A + V[a f ]fi c dt (54) 
+ y/^Hlde-**]^ dW{t) . 

and the expression for the homodyne current, 



i c {t) = e r (c, 



,-itj> 



Here we set 6 = r Im[/(— v 
atomic "jump operator" , 



+ y/eT~ m . (55) 
I{v)]/A and defined the 



G 



\/2vT 

n 2 



I{-v)a + I{+v)a)] 



(56) 



This operator as well as the laser heating and cooling 
rates, A± = roRe[/(±i/)]/2, depend on the function 
J(±i/), which is defined in Section [n] Eq. (fTtl)l . 

As in the previous section we consider a feedback force 
which is proportional to the measured signal, I c {t). Note 
that depending on values of 7(±z/), and the local oscilla- 
tor phase, 4>, the force is proportional to a linear combi- 
nation of (p) and (z) . For the derivation of the feedback 
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FIG. 6: a) Laser heating and cooling rates, A±, as a function 
of the laser detuning, A. b) Feedback heating and cooling 
rates, A^, for G = 1 and an optimized local oscillator phase, 
<j>. The parameters for the two plots are: e = 0.05, Q/T — 0.8, 
and v/Y = 0.1. 



master equation we follow the outline given in Section HV 
and obtain 



(i = — i(v+ 5)[tfa, n] + A-V[a](i, + A + I?[a t ]/i 

Z 06 



(57) 



Under the rotating wave approximation, which is valid 
for ro <C v, and by neglecting small shifts of the trap 
frequency, we end up with master equation @ given in 
Section in) 

Discussion: Fig. [S] shows the dependence of the four 
different rates A±, A{ b as a function of the detuning A. 
The cooling and heating rates which originate from the 
laser interaction, A± , correspond to the rates for ELC de- 
rived in Ref. |l8j . For the parameter regime > 1v and 
for blue detuning, A > 0, they lead to a minimal tem- 
perature for fi 2 ss \v(y — A). For red detuning, A < 0, 
the heating rate is larger than the cooling rate and with- 
out feedback the system does not reach a steady state. 
By adjusting the phase cj> the feedback loop always pro- 
vides additional damping, W — 

A f_b A fb 

> 0, which 

for A = and ^ <gC O is given by W = GTq. The noise 
added by the feedback loop r G 2 /8e, imposes a restric- 
tion on G if one is interested in low steady state energies. 
The combined effect of feedback and EIT laser cooling 
lead to a final temperature which is plotted in Fig. [3] in 
Section ELDl 



11 



VI. CONCLUSION 

In this paper we have shown that a continuous readout 
of the momentum of a single atom can be achieved by 
employing the high velocity sensitivity of the index of 
refraction in a driven A-system. The transparency effect 
for an atom at rest and the linear dependence of the index 
of refraction on the Doppler shift lead to a homodyne 
current linear in (p) and a minimal back action on the 
atomic motion, approaching the quantum limit for e — ► 1. 

By applying a force which is proportional to the mea- 
sured signal, feedback cooling for single ions can be real- 
ized. The cooling scheme is applicable in and outside the 
Lamb-Dicke regime with steady state temperatures well 
below the Doppler limit. From a fundamental point of 
view we want to point out, that in the proposed feedback 
scheme the measured current is fed back directly on the 
trap electrodes. Therefore, its implementation allows for 
a test of the theory of direct quantum feedback [| on an 
individual quantum system close to the ground state. 
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APPENDIX A: ADIABATIC ELIMINATION 

In this appendix we derive the stochastic master equa- 
tion (|38[l for the external density operator, fi c , and the 
expression for the signal I c (t) l|39|l . We start with the 
conditioned master equation given in Eq. (|28|l and apply 
the unitary transformation as defined in Eq. i|29|l . In the 
new basis the stochastic master equation can be written 
in the form 

dp c = (d + £„ + j) p c dt + ^/eT g H[<j ge e-^}p c dW(t) , 



(Al) 

where we divided the total Liouville operator into the 
three contributions, 

£i(p) = ~[H I ,p]-^(\e)(e\p + p\e)(e\) ) 

1 

£>v{p) = --[H E + H in t, p], 
n 

J(p) = rj g (a ge pa eg ) +TJ r (a re pa er ) , . 



The action of the recycling operators, Sf g>r , is defined in 
Section Hill Eq. (gTJ). 

In the following we write the total density operator in 
terms of the eigenbasis of Hi as, p c — J2i j Atj ® I*) (j\ 



with i,j E {+,—,£>}. Reinserting this decomposition 
into Eq. I|A1[) we obtain a set of coupled equations for 
the external operators jlij. By grouping the 9 elements 
pij into a single vector p, the resulting set of equations 
can be written in the form 



dp = (L/ + L„ + J)pdt- 



Tv{V-p}j pdW{t). 

(A2) 

The entries of the matrices L/, L„, J, S and the vector, 
V, can be derived in a straight forward (but lengthy) 
way by writing the operators in Eq. IjAlfl in terms of the 
states, |+), |— ) and \D). Note that the entries of L„ 
and J still contain superoperators acting on the external 
operators flij. 

The goal is to derive an effective equation for the 
population in the dark state, po := pdd, for the pa- 
rameter regime v, <C fl± . Formally we can com- 
bine both conditions and make a series expansion in 
the trap frequency, v. According to the structure 
of L/ we group the external operators p^j into the 
two vectors pi := {p D+ , fi+D, M-D-- P-d) T and p 2 := 

{p++, P--, Mh ,p h) T - By ordering the entries of p 

such that, p = (p2 , pf , po) T , the matrices Lj, J and 
L y have the block form 





From the structure of Lj and L y we see that for v — > 
we have pi ~ 0{v) and pi ~ 0{v 2 ). Therefore, up to 
second order in v the equation for po is 



dp = (L°Ji + K L p x + J^p 2 ) dt 



oi ,- 



7 02; 



+ ^/^~ g {v-pi - Tr{v- pi}p )dW(t) , 



(A3) 



with 



-jf (e** cos(0), cos(0), e** sin(6>), e"^ sin(0)) . 



The equations for p\ and pi are given by 
ft! = ^j^ + L^po+O^ 2 ), 



(A4) 



M2 



(Lj + J 2 ) p 2 + Li'p x + <D{v 



They can be integrated and up to the relevant orders of 
v we obtain the formal solution 



(^r^iVo + eV), 

(Lj + jy'Ll^Ljy^po + O^). 



mi 

fl2 



Resubstituting these expressions into Eq. I|A3() the result- 
ing equation can be written in the form 



dp = - i[h cS p - p hl s ]dt + X 2 TTZ(pp p)dt 



e ll±^H[pe-^]p dW(t), 



(A5) 
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with A = 2rjiygflL/n 3 , the non-hermitian operator 



^cff 



and the "recycling" term, 



K = 



TP 



1 



111 
TP 



(A6) 



(A7) 



In the formal expression for 1Z the inversion of operator 
is justified in the limit of a weak probe field, g <C £Il and 
r, r <C r. From the stochastic term in Eq. (|A5() we see 
that these conditions also maximize the signal strength 
for a given decoherence rate, A 2 T. 

In the original basis the evolution for [i c is given by 
the relation, 



dn c = Tn{Udjlo <g> \D)(D\U*}. 



(A8) 



Due to the overlap |(r|£>)| 2 = g 2 /fl 2 the action of U on 
the external states reduces to the action of the operator 
exp(ifjzg 2 /il 2 ). Therefore, the only effect of the basis 
transformation is the cancellation of the term Xgp/2 in 
the effective Hamiltonian, h c ff. 

For the expression of the measured signal I c (t) i|39|) we 
can simply repeat the calculations from above. Using the 
same notation as in Eq. (|A2|) it can be written as 



(A9) 



7 c (t) = eT g Tr{V-fi} 



We see that the first term already appeared in the 
stochastic master equation (|A2|) and can be evaluated 
along the same lines. By multiplying the resulting ex- 
pression by a factor A, we obtain the current given in 
Eq. JS3>- 



APPENDIX B: ADIABATIC ELIMINATION IN 
THE LAMB-DICKE LIMIT 

In this appendix we derive the conditioned evolution 
of the external atomic state in the Lamb-Dicke limit, 
VpiVl <C 1. We start with the stochastic Schrodingcr 
equation for the total wavefunction, |^), which includes 
the state of the atom as well as the state of the elec- 
tromagnetic environment. For the system Hamiltonian 
Hs (|52|) and the atom-field interaction described in Sec- 
tion ^0] it is given by 



dm 



-Ijf< ? -~|e)< e | j |*W/ 



(Bl) 



The noise increment operators, dB^(t), fulfill the Ito rules 
dBu^dB^^t) = 8(u — u')dt and correspond to the emis- 
sion of photons under an angle a — arccos(it) with re- 
spect to the z-axis. 



We decompose the total wave function in terms of the 
eigenstates of Ei as |*) = J2i |«), with i = +, — , D. 

For vanishing Lamb-Dicke parameters, Tjj — > 0, and af- 
ter some transient deviations the system evolves into the 
state, \^>) = \^d) <8> \D). The coupling between the exter- 
nal and internal degrees of freedom, -ffint, leads to finite 
contributions from the bright states which are of the or- 
der of the Lamb-Dicke parameter, |^±) ~ 0(fj). In the 
following we treat the coupling to the excited states in 
perturbation theory to derive an effective equation for 
\^ d) which is valid up to second order in fj. 

In the interaction picture with respect to the external 
Hamiltonian, He = hvcvd, the equation for the dark 
state wave function is 



■f m\*e)dt- 



Vf / duy/N(u)dBi(t)\V t 



(B2) 

with |* e ) = cos(6»)|* + ) + Bin(0)|*_). Note that in this 
equation we already used the assumption of a weak probe 
field, and set f2 = V^l + g 2 ~ Ql- Under the same 
assumption the equations for the bright states are given 
by 



d_ ( |* + ) 
dt { |*-> 



= M 



I*- 



where we defined the matrix 



M 



■cos 2 (9) 



sin(20) 



sin(20) 



\6) 



(B3) 



(B4) 



Up to first order in 77, the solution for the excited states 



1*+) 



2 



-M(t-s) 



z(s)ds 



cos(0) 
sin(0) 



I* 



D 



(B5) 

By inserting the time dependence of the position oper- 
ator, z{t) — (ae~ lvt + a) e lut ) / 'v2, we can evaluate this 
integral and obtain the evolution of the wavefunctions, 
\^±). This solution is then reinserted into Eq. (|B2|I to 
get the equation of the dark state wave function up to 
order fj 2 . As long as the final dynamics in the interac- 
tion picture is slow compared to the trap frequency, u, 
we can use the rotating wave approximation and neglect 
terms proportional to e ±l2vt . The resulting equation is 
then given by 



-2 ? 

V 9 



I(v)da) + I(-v)a)d \^ D )dt 



2 2 /*! 

' rg r / duy/N{uj C{t)dBl (t)|* D ) , 



(B6) 

where we set C(t) = (I(—v)a(t) + I(v)aJ(t)), and defined 
the function I(v) by 



I{v) := (cos(0),sin(0)) 



e- MT e~ iUT dT 



cos(0) 
sm(9) 



(B7) 
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Apart from the motional state of the atom, the wave- 
function I^d) still incudes the full state of the electro- 
magnetic environment. To obtain the conditioned dy- 
namics for the external density operator, /i c , we first de- 
compose the set of noise increment operators into the two 
contributions, 

dBl(t) = -= / duy/N&dBtit), 

dB\(t) = -JL= J e du^T)dBi(t) . 

As in Section ITTll the parameter, e, determines the frac- 
tion of the photons which are scattered into the mode of 
the probe beam. The increment operator dB* p (t) obeys 
the Ito rule, dBp^dB^t) = dt, and corresponds to the 
emission of photons, which are focused on the detector. 

A rigorous, but rather technical way to convert the 
stochastic Schrodinger equation into a stochastic mas- 
ter equation can be found in Ref. |26;, 30] . A convenient 
shortcut is, to first derive the unconditioned master equa- 
tion for /i = TrBji/{|4 f £))(4 f £)|} and then perform a par- 



tial 'unravelling' of this master equation as discussed in 
Ref. HI. For a phase <j> of the local oscillator we finally 
get the equation 

—1 9 

dfj, c = — i — - — Im[J(— v) + [a 1 a, ii c ]dt 

8 

—2 2 

+ (Re [/(-!/)]£> [a] +Re[7(+i/)]2?[o t ]) fi c dt 

+ \f^fn[c(t)e-^ c dw(t) 1 

(B8) 

The measured current has the form 



I c (t) = eY^f{C{t)e-^ + &(t)e«) c + ]f^fm ■ 

(B9) 

To make a comparison to results of Section IIIII and Sec- 
tion \W\ we introduce the rescaled function, 1(f), by set- 
ting I(v) = I(v)tt 4 /(8v 2 T). With this definition and 
a rescaling of the current we finally obtain Eqs. Ij54(l 
and |JSSJ. 
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